Collisional dynamics of ultracold OH molecules in an electrostatic 

field. 



Alexandr V. Avdeenkov and John L. Bohn 
JILA and Department of Physics, University of Colorado, Boulder, CO 

(February 2, 2008) 

Abstract 

Ultracold collisions of polar OH molecules are considered in the presence of 
an electrostatic field. The field exerts a strong influence on both elastic and 
state-changing inelastic collision rate constants, leading to clear experimental 
signatures that should help disentangle the theory of cold molecule collisions. 
Based on the collision rates we discuss the prospects for evaporative cooling 
of electrostatically trapped OH. We also find that the scattering properties 
at ultralow temperatures prove to be remarkably independent of the details 
of the short-range interaction, owing to avoided crossings in the long-range 
adiabatic potential curves. The behavior of the scattering rate constants is 
qualitatively understood in terms of a novel set of long-range states of the 
[OH] 2 dimer. 
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I. INTRODUCTION 



Polar molecules bring something entirely new to the field of ultracold physics. As com- 
pared to the neutral atoms that have been studied experimentally in the past, polar molecules 
possess extremely strong, anisotropic interactions. It has been speculated that dipolar inter- 
actions will lead to new properties in Bose-Einstein condensates |[T|-§[ or degenerate Fermi 
gases [Q. It has also been suggested that polar molecules in optical lattices may be useful in 
implementing quantum logic elements j^. On the experimental side, cold polar molecules 
may be produced in several ways, including photoassociation of two distinct alkali species 
PJ^, buffer-gas cooling or Stark slowing |T0|-[T2|. 

Regardless of the method or production, collisions of molecules are of paramount impor- 
tance in describing the properties of the gas. Collisions should also be interesting in their 
own right, as detailed probes of intermolecular interactions. Several features of the colli- 
sional dynamics of ground-state polar molecules, based on a simplified "toy" model, were 
discussed in . This model accounted for the interplay between dipole-dipole interactions, 
an external electric field, and states of different parities. The dipole-dipole interaction, which 
scales with intermolecular separation i? as renders cold molecule collisions completely 

different from cold atom collisions. This is because a l/i?^ interaction is characterized by 
energy-independent low energy cross sections in a// partial waves, not just in s- waves |T^-|TB[. 
The relatively strong, long-range interactions imply that molecules electrostatically trapped 
in weak-field-seeking states are generally susceptible to state-changing collisions that can 
rapidly deplete the trapped gas. The rates are in general far larger than those of magnetic 
dipolar transitions in stretched-state alkali atoms, owing largely to the relative strength of 



electric, as opposed to magnetic, dipolar interactions [13 



In this paper we address ultracold polar molecule collisions in a more realistic model, 
considering in detail the OH radical. This choice is motivated by the attractiveness of this 
molecule for Stark slowing from a supersonic jet [p!2| , p!7| . In particular, it has a ^11 ground 
state with a small A-doublet splitting, making it easily manipulated by modest-sized electric 
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fields. A full treatment of cold collisions is somewhat hindered by the fact that the OH-OH 
potential energy surface (PES) is poorly known, although it is known to be very deep and 
strongly anisotropic |I9|J20| . It is not even known, for example, whether OH molecules may 
suffer chemical reactions at ultralow temperatures. As a point of reference, it was recently 
suggested that the reaction F+H2 — HF + H may proceed at appreciable rates at ultralow 
temperatures, in spite of having a chemical barrier height of 700 K 

However, long-range dipole-dipole forces strongly dominate the scattering of OH 
molecules in their weak-field-seeking states. In this paper we will show that this arises 
from strong avoided crossings in the long-range adiabatic potential curves, which prevent 
the molecules from approaching close enough to one another for exchange potentials to 
become important. In this regard cooling and electrostatic trapping of OH molecules can 
provide a wealth of information on the long-range OH-OH interaction. Thus it appears pos- 
sible to understand a class of ultracold OH-OH collisions without detailed knowledge of the 
short range PES. This strategy would be an important stepping stone toward understanding 
the full problem of ultracold OH collisions. Strong- field seekers, by contrast, will in addi- 
tion experience the short-range interaction. The complete problem of exploring collisions 
of ultracold polar molecules might therefore most efficiently proceed by a two-step analysis, 
thus simplifying this very complicated problem. 

Accordingly we focus in this paper on the first step, namely, collisions of weak-field- 
seeking states. After some discussion of the relevant properties of OH molecules in Sec. II 
and their interactions in Sec. HI, we move on in Sec. IV to illustrate some prominent energy- 
and field-dependent features in elastic and inelastic cross sections. Mapping these features 
in experiments should help in unraveling the long-range part of this puzzle. We also present 
a simplified model of the long-range interaction, to help illustrate the basic physics behind 
the behavior of the cross sections. It will turn out that a new class of long-range bound 
states of the [OH] 2 dimer play a significant role in ultracold collisions of this molecule. 
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II. OH MOLECULE 



The OH molecule has a fairly complicated internal structure, incorporating rotation, 
parity, electronic spin, and nuclear spin degrees of freedom, which are further confounded 
in the presence of an electric field. We therefore begin by describing the structure of this 
molecule and the simplifications we impose to render our model tractable. 

Molecules cooled to sub-Kelvin temperatures by Stark slowing will be assumed to be in 
their electronic ^11 ground state, and v = vibrational ground state. In this state OH is 



an almost pure Hund's ) molecule, and has a dipole moment of 1.68 D ||19[. Spin- 

orbit coupling involving the lone electronic spin splits the ground state into ^H3/2 and ^Hi/2 
components, of which ^11^/2 is lower in energy and is therefore the state of greatest interest 
in ultracold collisions. In our model we take into account just the lowest rotational level 
of the corresponding ground state, J = 3/2. The energy of the first rotationally excited 
state with \J = 5/2, i7 = 3/2 > is about 84K higher in energy ||18| and we will neglect this 
and higher-lying states in our scattering calculations. Such states will, however, contribute 
rotational Feshbach resonances in realistic collisions. 

The isotopomer ^^O^H that we consider here has a nuclear spin of J = 1/2, which with 
a half- integer rotational quantum number defines the OH molecule as a boson. Thus we 
should take into account hyperfine structure to ensure the proper Bose symmetry. We will 
see below that the inclusion of hyperfine structure is also important in determining details 
of collision properties. The calculations in an electric field also require knowing the Stark 
splitting for OH molecules. Thus the Hamiltonian for the OH molecule in a field is 

H^^ = Hrot + Hfs + Hhfs + H field (1) 

Wave functions for the spatial degrees of freedom of the molecule are constructed in the 
usual way. Namely, in the zero-electric-field limit, eigenstates of parity e (= ±) are given 
by the Hund's case (a) representation: 

iJMjQe >= -^{iJMjQ > |AS > +e\JMj - > | ~ A - S >), (2) 
V 2 



where the rotational part is given by a Wigner function: 

,2J+ 1, 



UMrn >-- 



(3) 



and Q = |S + A| is the projection of the total electronic angular momentum on the molecular 
axis. The total spin of the molecule, F = J + I, with laboratory projection Mp, is then 
constructed by 



\{JI)FMFne >= |A > |5S > iJMjQe > \IMi >< FMf\JMjImj >, (4) 

Mj,Mi 

The matrix elements for the Hamiltonian (|l]) in this basis can be found elsewhere p3 



In compact form these matrix elements are 

< {JI)FQMF6\H''''\{J'l')F'Q'MF,e' >= (5n,3/25n',3/2^3/2,3/2 + 
<^C,l/2<^n',l/2-E'l/2,l/2 + ('^C,3/2<^C',l/2 + '^n,l/2'^n',3/2-£'3/2,l/2) X 



(5) 



X 



mj'][F][F']) 



VI N 1/2 



V 




In this expression yU is the molecular dipole moment, £ is the strength of the electric field, and 
Eq„q.' are matrix elements for the fine structure ifj-ot + -f^/s which can be found in p3| , p4| . 
These values depend on the rotational constant, spin-orbit coupling constant, hyperfine 
coupling constant, and A-doublet parameters of OH. All of these constants can be found 

HI. 



m 



Equation shows that, strictly speaking, the only good quantum number for the OH 
molecule is the projection of its angular momentum on the laboratory axis, Mp. However, 
for our present purposes it suffices to treat the quantum numbers as "almost good." For 
example, in view of the fact that OH is nearly a purely Hund's ) molecule, the 

coupling between Vt = 1/2 and f2 = 3/2 states is fairly weak. We account for this interaction 
perturbatively, by replacing the values E^i2 and -E1/2 by the eigenvalues of the 2x2 matrix 

\ 

(6) 



/ 



-£'3/2 -£'3/2,1/2 
y -£'3/2,1/2 Ei/2 y 
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keeping all other quantum numbers constant. 

Likewise, different values of the molecular spin J are mixed in a field, but this mixing is 
small in laboratory strength fields. The total spin F and the parity are far more strongly 
mixed. Accordingly, in practice we transform the molecular state to a field-dressed basis for 
performing scattering calculations: 

\{JI)FMFne;S >= ^ a{JFe)\{JI)FMFVte >, (7) 

JFe 

where a{JFe) stands for eigenfunctions of the Hamiltonian (|l]) determined numerically at 
each value of the field. We will continue to refer to molecular states by the quantum numbers 
J, F, and e, with the understanding that they are only approximately good in a field, and 
that is the appropriate molecular state. 

Figure 1 shows the Stark energies computed using all the ingredients described above. 
In zero field the energy levels are primarily determined by the A-doublet splitting between 
opposite parity states, whose value is A = 0.056cm~^. The alternative parity states, with 
e = —1 (f states) and e = +1 (e states) are shown in Fig. 1(a) and Fig. 1(b), respectively. 
These states are further split into hyperfine components with total spin F = 1 and F = 2. 
The Stark shift is quadratic for fields below the critical field So = A/2yU (^ lOOOV/cm for 
OH). For fields larger than Sq states with different parity are entirely mixed and the Stark 
effect transforms from quadratic to linear. In this case the molecular states are roughly 
equal linear combinations of the zero-field e = — and e = + states (compare Eqn. |^) p5| : 



\JMjT^ > Mj <0 
\JMjQ6 = ±l>={ (8) 

\JMj±n > Mj > 



III. OH-OH INTERACTION 

We will consider diatom-diatom scattering as two interacting rigid rotors in their ground 
rotational states. The complete Hamiltonian for the collision process can then be written 

H = Ti + T2 + //f ^ + H^"" + Vs + V^^ + V,, + Vmsj>, (9) 



where Tj and H^^ are the translational kinetic energy and internal motion of molecule 
i, including the electric field as in Eqn. ([^). Vg is the short- range exchange interaction, 
Vfj.fj. + Vqq + Vdisp are the dipole- dipole, quadrupole- quadrupole and dispersion long- range 
interactions respectively. Explicit expression for the dipole-dipole (oc l/i?/^) and quadrupole- 
quadrupole (oc interactions are given in Ref. Matrix elements for the the dipole- 



quadrupole interaction vanish for rigid rotor molecules in identical states [1^, hence will not 
be considered here. 

The anisotropic potential between two interacting rigid-rotor molecules is conveniently 



recast into a standard set of angular functions [p6| : 



Vs + Vf,^ + Vgq + Vdisp = V{uja, ujb, uj,R) = Y^ VaAa{uja, ^^b, uj), (10) 



where A = {La, Ka, Lb, Kb, L) and the angular functions are defined as: 

DitA,KA^A)D^^r,,K,MClM, (11) 



AA{u}A,i^B,^) = XI 

Ma,Mb,M 



^ La Lb L ^ 



yM^ Mb M j 



where uja,b = {0a,b, <Pa,b) are the polar angles of molecules A and B with respect to the 
lab-fixed quantization axis, and R = {R, u) is the vector between the center of mass of the 
molecules in the laboratory fixed coordinate frame. The indices Ka and Kb denote the 
dependence of the interaction on the orientation of the molecules about their own axes; in 
what follows we will ignore this dependence, setting Ka = Kb = 0. For the long-range 
part of the interaction this approximates the quadrupole moment of OH as cylindrically 
symmetric. 

The exchange potential Vs is very complicated, consisting of four singlet and four triplet 



surfaces pOj) is moreover poorly characterized. The most complete treatment of this 
surface to date computes the lowest-energy potential for each value of internuclear separation 
R |T^. This potential finds an extremely deep minimum at i? = 2.7 a.u. corresponding to 
chemically bound hydrogen peroxide, and a second, shallower minimum at = 6 a.u. due 
to hydrogen bonding forces. However, in cold collisions, the scattering cross sections are 
so sensitive to details of the short-range interaction that knowing the complete interaction 



probably would not help anyway. More importantly, as we will see below, collisions of the 
weak-field-seeking states are strongly dominated by the long-range dipole-dipole interaction. 
Therefore, we will use at small R simply the hydrogen-bonding part of the potential surface 
(see Fig. 13 of Ref. [0), and we will treat this part of the interaction as if it were isotropic. 
Finally, we will assert that the spin states of the OH molecules are in their stretched states, 
so that ordinary spin-exchange processes will not play a role in these collisions. 
We express the Hamiltonian in a basis of projection of total angular momentum, 

M = Mf, + Mf, + Ml, (12) 
Mf^=Mj^ + Mi^, (13) 

where Mp^, Mj-, Mj. are the projections of the full molecule spin, rotational motion and 
nuclear spin on the laboratory axis respectively for each molecule. Mi is the projection of 
the partial wave quantum number on the laboratory axis. In this basis the wave function 
for two molecules is described as: 

^■^ = ^ {|1 > ®|2 > ®\lMi >}^ X tfj^'^'^R), (14) 

1,2,/, Mj 

where {...}'^ is angular momentum part of this wave function and |i > is the wave function 
for each molecule as described by Eq.(|). 

Because the target and the projectile are identical bosons, we must take into account the 
symmetry of the wave function under exchange. The properly symmetrized wave function 
is then 

{|1 > ^12 > ^Mi ^y=i\^>^\^> ^1^^^ >} + i-^y i\^>^\^> ^1^^^^ >} (15) 

V2(l + S12) 

Using the expansion of the intermolecular potential (p!0|), the wave function (^2|), and 
taking into account the Wigner-Eckart theorem, we can present the reduced angular matrix 
element as 

< 12/Mi|||AA|||l'27'Mp >= (16) 
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f _^\LA+LB+Ji+Ji+J2+J2+M'p^ +M'p^ -n[-n'2+Mi-l (1 + gigi(-l)''-'^) (1 + g2g2(~l)^^) 

' ' 2 2 



X 



^ t/( La Ji \ J2 Lb J2 ^ ^ 



n[ 



\^ Mf, - M^, M^, - Mf> Ml - Ml, J 



V 



Vl'^ -VL2 
F, Fl\ ( I' 



y Mf, - Mf^ -Mf, Mf^ j 



La F, F[ 

y Mf^ - Mf[ -Mf^ Mf[ j 
L I \ ( I' L A 







Mi> Ml - Mv -Ml 

La -^1 -^1 I I Lb F2 F2 
I Jv Ji \ \ I J2' J2 

where / = 1/2 is each molecule's nuclear spin. 

The reduced matrix elements of the angular functions between symmetrized basis 
states (p!5D are 



< 12/M;||Ui|||l'27'Mi/ >= 



< 12/Mz|P^|||l'27'Mi/ > +(-1)' < 2llMi\\\A%\\\V2'l'Mv > 



;i + 5i,2)(l + 5i',2' 



(17) 



1 + f-l 



X- 



In practice, before each scattering calculation we transform the Hamiltonian matrix from this 
basis into the field-dressed basis defined by (J^). We solve the coupled- channel equations 



using a logarithmic derivative propagator method |22] to determine scattering matrices. 
Using these matrices we calculate total state-to-state cross sections and rate constants. 



IV. RESULTS 

This paper considers the scattering problem for OH molecules in an electrostatic field for 
cold and ultracold temperatures. We are interested in particular in the highest energy weak- 
field-seeking state of the ground rotational state, |(J, I)FMf, >= |(3/2, 1/2)22, 3/2, — >. 
This state is indicated by the heavy solid line in Figure 1. Since the quantum numbers J, 
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/, and n are the same for all the scattering processes we will consider, we will refer to this 
state by the shorthand notation \FMp,e) = |22, — >. 

The main novel feature of OH- OH scattering, as compared to atoms or nonpolar 
molecules, is the presence of the long-range dipole-dipole interaction and its dependence 
on the electrostatic field. Because these interactions strongly mix different partial waves, it 
is essential that we include more than one value of /. However, in the interest of emphasizing 
the basic underlying physics, we have included only the s- and d- partial waves. Sample 
calculations show that higher partial waves change the results only slightly at the energies 
considered. In this case, given the initial state with Mpi = Mf2 = 2, the only allowed 
values of the total projection are A4 = 2,3, 4, 5, 6. Among these channels only the one with 
A4 = 4 contains a contribution from s-wave scattering, and so will deserve special attention 
in what follows. In this case the total number of scattering channels for all allowed values 
oiM,is 208. 



A. Prospects for evaporative cooling 



One of the goals of the present work is to revisit the conclusions of Ref. concerning 
the effectiveness of evaporative cooling for electrostatically trapped molecules. To this end 
Figure 2 plots the elastic and state-changing inelastic rate constants versus field strength 
for two different collision energies, 100 fiK and IfiK. Here "elastic" refers to collisions that 
do not change the internal state of either molecule, while "inel" denotes those collisions 
in which one or both molecules are converted into any other states. These transitions are 
typically exothermic, leading to trap heating. Not all these collisions produce untrapped 
states, however. We find that the main contributions to the Kmei are given by processes 
in which quantum numbers F and/or Mp are changed by one. In particular, the process 
|22, — ) + |22, — ) |22, — > +|21,— > generally makes the largest contribution to Kinei, 
especially at high electric field. 

At low field the rates are nearly independent of field, but begin to evolve when the field 
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approximately exceeds the critical field Sq = A/2fi where the Stark effect changes from 
quadratic to linear. Above this field the rate constants exhibit oscillations as a function of 
field. These oscillations provide an experimentally variable signature of resonant collisions, 
meaning that mapping this field dependence should help in untangling details of the long- 
range OH-OH interaction. This is similar to the ability of magnetic-field Feshbach resonances 



in the alkali atoms to yield detailed scattering parameters [^^|. 

Following the example of ultracold atoms, we expect that evaporative cooling can proceed 
when the ratio of elastic to inelastic collisions, Kei/Kinei ^ 1- Figure 2 shows that this is 
hardly ever the case for large field values £ > except, perhaps, at very special field values 
where Ki^ei is at a minimum of its oscillation. Since the losses are dominated by exothermic 
processes, the ratio Kgi/Kinei in the threshold scattering limit scales as the ratio ki/kf of 
the incident and final channel wave numbers, as can be seen from the Born approximation. 
Thus at high electric fields, where the Stark splitting is large (hence kf is large), the ratio 
may become more favorable. In our calculations this apparently happens for fields above 
10^ V/cm. 

For fields below So ~ 1000 V/cm, a favorable ratio of Kei/Kinei is only somewhat more 
likely. For fields this low, however, the maximum depth of an electrostatic trap is ~ 8mK, as 
given by the magnitude of the Stark shift (Fig. 1). The temperature of the trapped gas must 
therefore be well below this temperature. In the example of a 100 fiK gas, (Fig. 2a), the ratio 
Kei/Kinei Hiay indeed be favorable. However, if the gas is cooled further, say to 1 fiK (Fig. 
2b), this ratio becomes less favorable again. This is because of the Wigner threshold laws: 
the exothermic rate Kinei is energy-independent at low energy, while the elastic scattering 
rate plummets to zero as the square root of collision energy. Thus, in general, evaporative 
cooling seems to be viable only over an extremely limited range of temperature and field 
range for the OH molecule, if at all. We therefore reiterate the message of Ref. |T3|, and 



recommend that cold OH molecules be trapped by a far-off-resonance optical dipole trap, 
in their lowest energy \F\Mp\,e) = |11,+) states. 

At this point we emphasize an essential difference between evaporative cooling of electro- 
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statically trapped polar molecules and of magnetically trapped paramagnetic molecules. For 
polar molecules the transition from weak- to strong-field seeking states is always exothermic, 
even in zero applied field. This is because the lower member of a A-doublet is always strong- 
field-seeking (e.g. Fig. 1). For paramagnetic molecules, by contrast, weak- and strong-field 
seeking states can be nearly degenerate at low magnetic field values (e.g., ^"^02 as discussed 
in Refs. [|^,^. In the present case, OH is also paramagnetic and hence could in principle 
be magnetically trapped. For example, the low-energy states with {FMpe) = |11, +) might 
be suitable candidates. The infiuence of electric dipole interactions on evaporative cooling 
of magnetically trapped OH has yet to be explored. 

B. Analysis of the long-range interaction 

The general behavior of the rate constants in Fig. 2 can be explained qualitatively by 
simplifying our model even further, to a case that contains only the essential ingredients: 
the dipole-dipole interaction, the A-doublet, and an electric field |T3|. Roughly speaking, the 
electric field has two effects on the molecules: 1) it mixes molecular states of opposite parity, 
thus creating induced dipole moments; and 2) the resulting dipole-dipole interaction strongly 
couples scattering channels with different partial waves, leading to long-range couplings 
between two molecules. 

As a starting point in this analysis. Figure 3 breaks down the elastic and inelastic rates 
into their contributions from different values of the total projection of angular momentum 
Ai. This is done for the rates calculated at an energy E = 100 fiK, from Figure 2(a). 
In both elastic and inelastic scattering, the rates are dominated by the contribution from 
A4 = 4, which, it will be recalled, is the only value of A4 that incorporates s partial waves 
in the present model. We will accordingly consider only this case in what follows. 

The model used to obtain the results in Figs 2 and 3 consists of 32 channels for the 
block of the Hamiltonian matrix with Ai = 4. To simplify the analysis of this block even 
further, we focus on the sub-Hamiltonian with fixed quantum numbers F = Mp = 2 for 
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each molecule. This reduces the effective Hamiltonian to six channels: there are three non- 
degenerate thresholds E,,-^^.^ corresponding to different possible values of the field-dressed 
parity quantum number Si = ±. For each of these three thresholds there are two channels, 
corresponding at large R to s- and d- partial waves. 

The simplified six-channel Hamiltonian matrix then consists of 3 x 3 blocks y"' param- 
eterized by partial wave quantum numbers /, /': 



H 



^ yOO T>02 \ 

diag 



(18) 



Here the diagonal components Vf-ag include the parity thresholds and the centrifugal inter- 
actions, 



E_ 



2mB? 



V 


















where the electric-field-dependent thresholds are given by E^-^^^^ 
{si -\- £2)AvT+~P^/2, in terms of the dimensionless parameter 



(19) 



E+ + E_ - 



(20) 



that relates the electric field strength £ to the zero-field lambda-doublet splitting A = 
— The simplified field-dependent dipole-dipole interaction term V^'' is readily pa- 
rameterized in the field-dressed basis as 



II' 



e V2k 1 

V2k -^/2k 
1 -V2k k^ 



;i + F)i?3' 



(21) 



whose coefficient C^'', which is independent of both R and the electric field, is given by 



' I' 2 I \ Q?Ml{J{J + 1) + F{F + 1) - /(/ + 1)) 







2(J(J+1)F(F + 1))2 



(22) 
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Notice that the dipole-dipole interaction vanishes for s- waves, so that V^''" = 0. 

Within this simphfied model we will refer to the scattering channels by the parities ei and 
62 of the two molecules, along with the partial wave quantum number /. Thus the incident 

channel for weak-field-seekers will be denoted \eie2,l) = \ ,0). Recall that all other 

quantum numbers (J, /, Q, F, Mp) are assumed to have fixed values for each molecule. 

The explicit field-dependence in the coupling matrix (pi]) explains qualitatively the be- 
havior of our ultracold weak-field-seeking molecules, which have incident quantum number 
e = —. For zero electric field {k = 0), there is no direct dipole-dipole coupling between 
identical molecules. There is, however, an off-diagonal coupling to different channels with 
opposite parity, as can be seen in the form of the Hamiltonian (|l^). This interaction brings in 
the dipole-dipole coupling in second order, contributing an effective dispersion-like potential 
Cq^^ /R^, with a coefficient 

OS" - ^ « ^ (23) 

for both s- and d- partial waves. For the OH molecule this effective coefficient is ~ 4 x 10^ 
atomic units, far larger than for the alkali atoms that are familiarly trapped. Thus, even in 
zero external field the effective interaction strength of polar molecules is quite large. This 
may imply the breakdown of the contact-potential approximation to describing Bose- Einstein 
condensates of polar molecules, even when their dipoles are not aligned by an external 
field We note that the quadrupole-quadrupole interaction is relatively unimportant, 

becoming larger than this effective dispersion interaction only when i? >^ 3 x 10^ a.u. 
When the field is switched on, the s-wave channels undergo a qualitative change. Now 

the incident channel | , 0) sees a direct coupling to its d-wave counterpart | ,2), via 

the matrix element V^'=0'''=2 (x This perturbation generates an far stronger effective 

long-range potential ot the form C'/'/R', with 

where 1 = 2. Thus the electric field is able to completely alter the character of the inter- 
molecular interaction. 
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For d-wave collisions, the dipole-dipole coupling is direct, but not in the limit of zero field, 
where the molecules are in parity eigenstates. At low fields (A; <^ 1, where the Stark effect 
is quadratic), the diagonal coupling V^^ oc + k"^), is small. In this limit the molecules 

are nearly in parity eigenstates, hence do not "know" that they have dipole moments. At 
larger fields this interaction grows in scale, thus "activating" the dipoles. This is why the 
rate constants shown in Figs 2,3 begin to evolve at fields near S = Sq. It is also why the 
contributions from all angular momentum projections except the one with A4 = contribute 
only weakly to scattering at low field. The channels with Ai ^ 4 are all of d-wave character, 
hence obey a threshold law cr oc i?^ at low fields, then evolve to a a oc const threshold law 
at larger fields. 



C. Large-field oscillations and long-range states of the [OH] 2 dimer 

At fields larger than the critical field Sq, the rate constants in Fig. 2 exhibit oscillations 
with field. Significantly, these occur only when the projection of total angular momentum 
Ai = A, which is the only case in which s- and d- partial waves are mixed (Fig. 3). To 
understand this oscillating behavior of cross sections we show in Fig. 4(a) the adiabatic 
potential curves in the simplified 6-channel model (^). In the case shown the field is £^ = 10"^ 
V/cm. In this figure a strong avoided crossing can be seen at i? ~ 60 a.u., corresponding 

to the crossing of the attractive |£:i£:2, = I 5 0) channel with the repulsive | h, 2) 

channel. The strong dipole-dipole interaction between these different partial waves creates 
the adiabatic potential shown as a heavy black line and labeled Uq. 

This potential curve supports bound states of the [OH] 2 dimer. These bound states 



are of purely long-range character, similar to the long-range states of the alkali dimers |31 



that have been used in photoassociation spectroscopy studies of ultracold collisions p2| , |33 
Moreover, in the case of the [0II]2 states the shape of the potential Uo, hence the energies of 
the bound states, are strongly subject to the strength of the applied electric field. The curve 
in Fig. 4(a) in fact possesses no bound states in zero field, but five by the time the field 
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reaches 10^ V/cm. More realistic adiabatic potentials are of course more elaborate, as shown 
in Fig. 4(b) for the more complete Hamiltonian that includes hyperfine levels. Nevertheless, 
in this figure, too, can be seen adiabatic potential wells that will support bound states. 

The significance of these curves is twofold: the crossing is very adiabatic, implying that 
coupling to lower-energy channels is weak, and that therefore the cross sections depend only 
weakly on details of the short-range potentials. This we have indeed verified by altering the 
short-range potential in the full calculation. 

Additionally, as the field strength grows and the potential becomes deeper, new bound 
states are added to the potential, causing scattering resonances to appear. This is the 
cause of the oscillations observed in the rate constants in Fig. 2. To illustrate this, we 
reproduce in Fig. 5 the complete elastic scattering rate constant (solid line), along with the 
same quantity as computed in the simple six-channel model (dashed line). The qualitative 
behavior is nearly the same, namely, oscillations appear at fields above Sq. Moreover, the 
arrows in the figure indicate the values of the field for which a bound state of Uq coincides 
with the scattering threshold. These fields correspond fairly well to the peaks, although they 
are somewhat offset by coupling to other channels. Nevertheless, this simple picture clearly 
identifies the origin of these oscillations with the existence of long-range bound states. 

These resonant states are not Feshbach resonances, since there is no excitation of internal 
states of the molecules; nor are they shape resonances in the usual sense, since there is no 
barrier through which the wave function must tunnel. Instead, they are the direct result of 
altering the interaction potential to place a bound state exactly at threshold [|14|. Probing 
these states through direct scattering of weak-field-seeking states should reveal details about 
the long-range OH-OH interaction, making possible a comparison of theory and experiment 
without the need to fully understand the short-range [OH] 2 potential energy surface. 
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V. CONCLUSION 



In this paper we theoretically investigated ultracold collisions of ground state polar di- 
atomic molecules in an electrostatic field, taking OH molecules as a prototype. Focusing on 
weak-field-seeking states, we have strengthened the suppositions in Ref. [0 that long-range 
dipolar interactions drive inelastic scattering processes that are generally unfavorable for 
evaporative cooling of this species. However, at electric fields above a characteristic value 
So, oscillations occur in both elastic and inelastic collision rates, implying that a regime may 
be found where the ratio Keiasud KineiasUc is favorable for cooling. Even though evaporative 
cooling may be difficult, the inelastic rates may nevertheless prove useful diagnostic tools 
for cold collisions of these molecules. The Stark slowing technique provides a means of 
launching a bunch of molecules toward a stationary trapped target, i.e., of making a real 



scattering experiment |]T2|,|T7|. 

For actual trapping and cooling purposes, for instance as a means of producing molecular 
Bose-Einstein condensates or degenerate Fermi gases, it seems likely that the molecules must 
be trapped in their strong-field seeking states. Collisions of these species will present their 
own difficulties, since they will depend strongly on the short-range part of the potential 
energy surface. However, the scattering length for OH-OH scattering may be determined 
by photoassociation spectroscopy to the long-range bound states we have described above. 



This will be analogous to the determination of alkali scattering lengths |p2| , |33| , except that 
microwave, rather than optical, photons will be used. The detailed properties of the long- 
range [OH] 2 states, and prospects for using them in this way, therefore deserve further 
attention. 

This work was supported by the National Science Foundation. We acknowledge useful 
discussions with J. Hutson and G. Shlyapnikov. 
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FIG. 1. The Stark effect in ground state OH molecules, taking into account hyperfine splitting, 
(a) shows the states that have odd parity e = — in zero electric field (f states), whereas (b) shows 
those of even parity (e states). The weak-field-seeking state with quantum numbers F = Mp = 2, 
the subject of this paper, is indicated by the heavy solid line. Note that states with Mp = ±|Mi?| 
are degenerate in an electric field. 
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FIG. 2. Rate constants versus electric field for OH-OH collisions -with molecules initially in 
their \FMp,£) = |22— > state. Shown are the collision energies WOfiK (a) and lfj,K{h). Solid lines 
denote elastic scattering rates, -while dashed lines denote rates for inelastic collisions, in -which one 
or both molecules changes its internal state. These rate constants exhibit characteristic oscillations 
in field -when the field exceeds a critical field Sq ~ 1000 V/cm. 
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FIG. 3. Elastic (a) and inelastic (b) rate constants versus electric field for the same circum- 
stances as in Figure 2(a). The rates are separated into contributions from different values of A4, 
the projection of total angular momentum on the laboratory z-axis. 
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FIG. 4. Adiabatic potential energy curves. The curves in (a) correspond to the simplified 
six-channel model described in the text, and show a long-range potential well (labeled Uq) that 
can hold bound states of the [0H]2 dimer. The curves in (b) are those for the more complete 
calculation that includes hyperfine structure. 
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FIG. 5. Elastic rate constants versus field, as in Figure 2. The solid line reproduces the elastic 
rate constant from figure 2. The dashed line is an approximate elastic rate constant based on 
the simplified six-channel model described in the text. The arrows indicate values of the electric 
field at which bound states of the long-range potential Uq (Figure 4) coincide with the scattering 
threshold. 
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